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. Abstract 

' In this paper we present a constructive method to control the bihnear Schrodinger equation via two 

controls. The method is based on adiabatic techniques and works if the spectrum of the Hamiltonian 
u: admits eigenvalue intersections, and if the latter are conical (as it happens generically). We provide 

sharp estimates of the relation between the error and the controllability time. 
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I. Introduction 

In this paper we are concerned with the problem of controlling the Schrodinger equation 

> 

S; t^=lHo + f2u,{t)Hk]tP{t). (1) 

cn : V fc=i / 

^ ■ Here ^ belongs to the Hilbert sphere S of a complex separable Hilbert space 7i and Hq, . . . , Hm 
^ are self-adjoint operators on T-L. The controls ui,...,Um are scalar- valued and represent the 

action of external fields. Hq describes the "internal" dynamics of the system, while Hi, ... , Hm 
■ the interrelation between the system and the controls. 
c3 ! The reference model is the one in which Hq = —A + Vo(x), Hi = Vi{x), where x belongs 

to a domain D C M" and Vq, . . ■ ,Vm are real functions (identified with the corresponding 
multiplicative operators). However, equation ([D can be used to describe more general controlled 
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dynamics. For instance, a quantum particle on a Riemannian manifold subject to extemal fields 
(in this case A is the Laplace-Beltrami operator) or a two-level ion trapped in a harmonic 
potential (the so-called Eberly and Law model [|2l, [|71). In the last case, as in many other 
relevant physical situations, the operator Hq cannot be written as the sum of a Laplacian plus a 
potential. 

The controllability problem consists in establishing whether, for every pair of states ipo and 
^pi, there exist controls Uk{-) and a time T such that the solution of ([U) with initial condition 
ipiO) = ipo satisfies ^{T) = ipi. The answer to this question is negative when "H is infinite- 
dimensional. Indeed, Ball, Marsden and Slemrod proved in [O a result which implies (see [301) 
that equation ^ is not controllable in (the Hilbert sphere of) 1-L. Hence one has to look for 
weaker controllability properties as, for instance, approximate controllability or controllability 
between the eigenstates of Hq (which are the most relevant physical states). However, in certain 
cases one can describe quite precisely the set of states which can be connected by admissible 
paths (see El, [El, [|23l). 

In [fT2l| an approximate controllability result for ([T) was proved via finite-dimensional geomet- 
ric control techniques applied to the Galerkin approximations. The main hypothesis is that the 
spectrum of is discrete and without rational resonances, which means that the gaps between 
the eigenvalues of should be Q-linearly independent. Another crucial hypothesis is that the 
operator Hi couples all eigenvectors of Hq. This result has been improved in []9l where the 
hypothesis of Q-linear independence was weakened. Similar results have been obtained, with 
different techniques, in [|22l (see also fl20l. [l23l). 

The practical application of the results discussed above entails three main difficulties: 

• In most cases the techniques used to get controllability results do not permit to obtain (even 
numerically) the controls necessary to steer the system between two given states. 

• Even in the cases in which one can get the controls as a byproduct of the controllability result, 
they happen to be highly oscillating and hence they can be difficult to implement, depending on 
the experimental conditions. Roughly speaking, since one should move in an infinite dimensional 
space with only one control, one should generate many iterated Lie brackets. This is particularly 
evident in the papers [|9l, [[T2l . where the use of Galerkin approximations permits to highlight 
the Lie algebra structure. 

• Explicit expressions of time estimates, for the norm of controls and for their total variations 
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are extremely difficult to obtain. (For some lower bounds of controllability time and estimates 
of the norm of the controls see [|9l.) 

In most of the results in the literature only the case m = 1 is considered. In this paper we 
study the case m = 2 and we get both controllability results and explicit expressions of the 
external fields realizing the transition. The system under consideration is then 

with H{ui, U2) = H0+U1H1+U2H2. The idea is to use two slowly varying controls and climb the 
energy levels through conical intersections, if they are present. Conical eigenvalue intersections 
have been used to get population transfers in the finite dimensional case in IfTOl . |fT3l . [[T9l . [|29ll . 
||3TI . Some preliminary ideas given in the present paper can be found in |[2l, where a specific 
example (which is a version of the Eberly and Law model) is analyzed, and in ifTTI . The main 
ingredients of our approach are the following: 

• The adiabatic theorem that, in its rougher form, states the following: let \{ui, U2) be an eigen- 
value of H{ui,U2) depending continuously on (mi,M2) and assume that, for every mi,M2 G K 
(K compact subset of M^), \{ui,U2) is simple. Let 0(^1,^2) be the corresponding eigenvector 
(defined up to a phase). Consider a path (mi,M2) : [0,1] — )■ K and its reparametrization 
{ul{t) , u^it)) = {ui{et),U2{et)), defined on [0, l/e]. Then the solution tpeit) of the equation 
■d^ = (ffp + u\{t)Hi + ul{t)H2)i!e{t) with initial condition V'e(O) = 0(ui(O), ^2(0)) satisfies 

\\i,,{l/e) - e'U{u\{l/e) ,ul{l/e))\\<Ce (2) 

for some d = '13(e) G M. This means that, if the controls are slow enough, then, up to phases, the 
state of the system follows the evolution of the eigenstates of the time-dependent Hamiltonian. 
The constant C depends on the gap between the eigenvalue A and the other eigenvalues. 

• The crossing of conical intersections. Generalizations of the adiabatic theory guarantee that, 
if the path (ni(-), U2(-)) passes (once) through a conical intersection between the eigenvalues 
Ao < Ai, then 

||^,(l/e) -e^%K(l/e),u|(l/e))|| < CV^ (3) 

where i^eit) is the solution of the equation = (Hq + u\{t)Hi + u|(t)if2)^e(^) with initial 
condition (/'^(O) = 0o(wi(O), m|(0)) and tpQ, 0i are the eigenvectors corresponding respectively 
to the eigenvalues Aq, Ai (see [TT^). Figure [H illustrates a closed slow path in the space of 
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Fig. 1. A slow path climbing the spectrum of H{ui,U2), plotted in function of (ui,U2)- 

controls producing a transition from the eigenvector corresponding to the eigenvalue Aq to the 
eigenvector corresponding to A2 by crossing two conical singularities. Notice that the path could 
not be closed if only one control was present. Indeed u(t) would pass back and forth through the 
same singularity and the trajectory would come back to the original state. One of our main results 
is that choosing special curves that pass through the conical singularity the estimate in (|3]) can be 
improved by replacing y/e by e. Hence if some energy levels Aq, . . . , \k of the spectrum of H are 
connected by conical singularities, then one can steer, in time \/e, an eigenstate corresponding 
to Aq to an eigenstate corresponding to Afc with an error of order e. 

• The behavior of the eigenstates in a neighborhood of conical singularities. If the path 
(mi(-), M2(-)) is a piecewise smooth curve with a vertex (i.e. discontinuity at the level) at 
the conical singularity, the state of the system evolves with continuity, while the eigenstates 
corresponding to the degenerate eigenvalue are subject to an instantaneous rotation. The angle 
made by the path at the vertex can be used to control the splitting of probabilities between the 
two energy levels (see Figure d). This splitting phenomenon has already been described and 
exploited for controllability purposes on a two-dimensional system in [fTTI . 

The ideas introduced above lead to the following result: if the energy levels Aq, . . . , Am are 
connected by conical singularities, then the system is approximately spread controllable, i.e., 
for every given e > and po, . . . ,Prn > such that J^^oPf ~ there exists a control u 
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defined on [0, k + 1 phases 'do, . . . , -(9^ G M, and a trajectory corresponding to u satisfying 
■0(0) = 00 and ||'0(l/£:) — Yl^j=iPj^^'^^^ji'^^)\\ ^ ^- Moreover the control can be taken of 
the form u(t) = 'j{et), where 7 : [0, 1] is characterized explicitly. Hence the method 

provides precise time estimates in relation with the required precision. The method cannot be 
easily reversed, in order to explicitly characterize paths steering a state which is spread on several 
eigenstates to a single one. The difficulty lies on the loss of information about the relative phases 
during adiabatic evolution. 

We finally remark that systems for which the method can be applied are rather frequent. 
Indeed intersections of eigenvalues are generically conical for Hamiltonians of the form —A + 
Vo + uiVi + U2V2, as explained in Section HIl 

The structure of the paper is the following. In Section III we introduce the framework and we 
state the main result. In Section UlI] we recall the time adiabatic theorem and some results on the 
regularity of eigenvalues and eigenstates of parameter-dependent Hamiltonians. In Section HVl we 
deepen our analysis of conical intersection; in particular, we state and prove a sufficient condition 
for an intersection to be conical. Section |V] is devoted to the construction of some special curves 
along which we can obtain our controllability result, while the proof of the main theorem is the 
subject of Section |VIl Finally, in Section IVIIl we show that the same controllability result holds 
also for more general curves than those presented in Section IVl 



January 12, 2013 



DRAFT 



6 



II. Definitions and main results 
We consider the Hamiltonian 

H{u) = Ho + uiHi + U2H2, 

for u = (^1,^2) G M^. From now on we assume that H(-) satisfies the following assumption: 

(HO) Hq is a self-adjoint operator on a separable Hilbert space H, and Hi and H2 are bounded 
self-adjoint operators on "H. 

When necessary, we also make the following assumption on the Hamiltonian H(-): 

(HI) There exists an orthonormal basis {Xj}j of the Hilbert space such that the matrix 

elements {Xj,HoXk), {Xj,HiXk) and {xj,H2Xk) are real for any j,k. 

Remark 2.1: Hypothesis (HI) ensures that, with each u and each eigenvalue of H(u) (counted 
according to their multiplicity), it is possible to associate an eigenstate whose components with 
respect to the basis {Xj}j are all real. 

A typical case for which (HO) and (HI) are satisfied is when Hq = —A + V, where A is the 
Laplacian on a bounded domain C M'' with Dirichlet boundary conditions, V E L°°(r2,R), 

= L^(f2,C), and Hi,H2 are two bounded multiplication operators by real valued functions. 
In this case the spectrum of Hq is discrete. 

The dynamics are described by the time-dependent Schrodinger equation 

= ff(u(t))^(t). (4) 

Such equation has classical solutions under hypothesis (HO), u(-) piecewise and with an 
initial condition in the domain of Hq (see ll26ll and also f3|). 

We are interested in controlling @ inside some portion of the discrete spectrum of H(u). 
Since we use adiabatic techniques, the structure of the spectrum shall satisfy some particular 
features: roughly, the portion of discrete spectrum we consider must be well separated from its 
complement in the spectrum of the Hamiltonian, and this property must hold uniformly for u 
belonging to some domain in R^. 

All these properties are formalized by the following notion: 

Definition 2.2: Let a; be a domain in M^. A map S defined on to that associates with each 
u G w a subset S(u) of the discrete spectrum of H(u) is said to be a separated discrete spectrum 
on CO if there exist two continuous functions /i, /2 : a; — )■ M such that: 
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• /i(u) < /2(u) and S(u) C [/i(u), /2(u)] Vu G u. 

• there exists F > such that 

inf inf dist(A,[A(u),/2(u)]))>r. 

uGo; AeSpcc(_f/(u))\E{u) 

Notation From now on we label the eigenvalues belonging to S(u) in such a way that we 
can write S(u) = {Ao(u), . . . , Afc(u)}, where Ao(u) < ■■■ < Afc(u) are counted according to 
their multiplicity (note that the separation of S from the rest of the spectrum guarantees that 
k is constant). Moreover we denote by 0o(u), . . . ,0fc(u) an orthonormal family of eigenstates 
corresponding to Ao(u), . . . , Afc(u). Notice that in this notation Aq needs not being the ground 
state of the system. 



Definition 2.3: Let S be a separated discrete spectrum on tu. We say that dH) is approximately 
spread-controllable on S if for every u", E cu such that S(u°) and S(u^) are non-degenerate, 
for every e {(f)o{u^), . . . , 0fc(u°)}, p G [0, l]'^^-'^ such that Yl'i=oP'i = 1' ^^'^ every e > there 
exist T > 0, ^9o, • • • , ^^fc e M and a piecewise control u(-) : [0, T] -)> such that 

k 

||^(T)-^p,e^''^-0,(ui)|| <£, (5) 

j=0 

where is the solution of & with ^'{0) = (p. 

Our techniques rely on the existence of conical intersections between the eigenvalues. Conical 
intersections constitute a well-known notion in molecular physics. They have an important role 
in the Bom-Oppenheimer approximations (see for instance [[H, ifTSl . [|27ll . where they appear 
for finite dimensional operators). In the finite dimensional case they have been classified by 
Hagedom [[T4|. 

A unified characterization of conical intersections seems to be missing. The following defini- 
tion meets all the features commonly attributed to them. 

Definition 2.4: Let H{-) satisfy hypothesis (HO). We say that u G is a conical intersection 
between the eigenvalues Xj and A^+i if Aj(u) = Aj+i(u) has multiplicity two and there exists a 
constant c > such that for any unit vector v G and t > small enough we have that 

Aj+i(u + tv) - Aj(u + tv) > ct . (6) 
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It is worth noticing that conical intersections are not pathological phenomena. On the contrary, 
they happen to be generic in the following sense. Consider the reference case where T-L = 

L\nx), Ho = -a + Vq: d{Ho) = H\nx) r^ H^i^x) L'^i^X), Hi = Vi, H2 = V2, 

with n a bounded domain of R"^ for some e N and Vj G C°(1],M) for j = 0, 1,2. Then, 
generically with respect to the pair (^i, V2) in C°{Q, R) x C°(fi, M) (that is, for all pairs, {Vi, V2) 
in a countable intersection of open and dense subsets of C°(r2, M) x C°(r2, M)), for each u G 
and A G M such that A is a multiple eigenvalue of H0 + U1H1 + U2H2, the eigenvalue intersection 
u is conical. 

In order to check that this is true, we can apply the transversal density theorem (see [[I] 
Theorem 19.1]) with ^ = C%n,R) x C°(fi,M), X = R"^, Y = C%n,R), p{{VuV2),u) = 
Vq + uiVi + U2V2, and 

W = {V e C°(fi, M) I -A + 1/ : H\n, C) n H^{n, C) ^ L'^{n, C) has multiple eigenvalues}. 

The covering of W by manifolds of codimension two is obtained in ll28l . based on the properties 
proved in [|6l (see also ifTTl ). We obtain that, generically with respect to {Vi,V2), the intersection 
of ^((^1, V2), M^) with W is transverse. Equivalently said, generically with respect to (Vi,V2), 
for every u G and A G M such that A is a multiple eigenvalue of —A + Vq + uiVi + U2V2, 
for every {vi, V2) G \ {0}, the line {{ui + tvi)Vi + {u2 + tv2)V2 | t G M} is not tangent to 
W, i.e., the eigenvalue intersection u is conical. 

Moreover, each conical intersection (mi,M2) is structurally stable, in the sense that small 
perturbations of Vq, Vi and V2 give rise, in a neighborhood of u, to conical intersections 
for the perturbed H. Structural stability properties can be proved without resorting to abstract 
transversality theory, as will be shown in Section |V] Theorem 15.101 

Our main result is the following: it states that spread controllability holds for a class of systems 
having pairwise conical intersections, providing in addition an estimate of the controllability time. 
As a byproduct of the proof, we will also get an explicit characterization of the motion planning 
strategy (the path 7(-) below). 

Theorem 2.5: Let H{u) = Hq + uiHi + U2H2 satisfy hypotheses (HO)-(Hl). Let S : u h-^ 
{Ao(u), . . . , Afc(u)} be a separated discrete spectrum on u; C and assume that there exist 
conical intersections Uj E co, j = 0, . . . , k — 1, between the eigenvalues Aj, Xj+i, with Xi{uj) 
simple if / ^ j,j + 1. Then, for every u° and such that S(u°) and S(u^) are non-degenerate. 
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for every e {0o(u°), . . . , 0fc(u°)}, and p e [0, 1]''+^ such that Yl'i=oP'i = 1' ^^^^^ C > 
and a continuous control 7(-) : [0, 1] — with 7(0) = u° and 7(1) = u^, such that for every 

e>0 

k 

11^(1/^) 

j=0 

where tp{-) is the solution of © with ip{0) = 0, u(t) = 7(et), and "i^o, • • • , ''^fc ^ are some 
phases depending on e and 7. In particular, dH) is approximately spread controllable on S. 

III. Survey of basic results 

A. The adiabatic theorem 

One of the main tools used in this paper is the adiabatic theorem ([[8]|, ifTSl . [|2T]|. [[24|): here we 
recall its formulation, adapting it to our framework. For a general overview see the monograph 
ll27l . We remark that we refer here exclusively to the time-adiabatic theorem. 

The adiabatic theorem deals with quantum systems governed by Hamiltonians that explicitly 
depend on time, but whose dependence is slow. While in quantum systems driven by time- 
independent Hamiltonians the evolution preserves the occupation probabilities of the energy 
levels, this is in general not true for time-dependent Hamiltonians. The adiabatic theorem states 
that if the time-dependence is slow, then the occupation probability of the energy levels, which 
also evolve in time, is approximately conserved by the evolution. 

More precisely, consider h{t) = Hq + ui{t)Hi + U2{t)H2, t E I = [to,tf], satisfying (HO), 
and assume that the map t {ui{t) , U2{t)) belongs to Assume moreover that there exists 

C such that (ui{t) , U2{t)) E u for all t G / and S is a separated discrete spectrum on to. 

We introduce a small parameter e > that controls the time scale, and consider the slow 
Hamiltonian h{et), t E [to/e^tf/e]. The time evolution (from to/e to t) U^{t,tQ/e) generated by 
h{e-) satisfies the equation i-^U^{t,to/e) = h{et)U'^{t,to/e). Let r = et belong to [to,tf] and 
To = to; the time evolution f/^(r, tq) := U^{t / e.tq/ e) satisfies the equation 

ie^U'{r,r,) = h{r)W{r,ro). (7) 

Notice that t/^(r, tq) does not preserve the probability of occupations: in fact, if we denote by 
P*(r) the spectral projection of /i(r) on Il(u(r)), then P^{t)U^{t,To) is in general different 
from f/^(r, ro)P*(ro). 
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Let us consider tlie adiabatic Hamiltonian associated with S: 

ha{r) = h{r) - ieP,{r)P,{r) - ieP^{T)P^{T), 

where P^{t) = id — P*(r) and id denotes the identity on "H. Here and in the following the time- 
derivatives shall be intended with respect to the reparametrized time r. The adiabatic propagator 
associated with ha{T), denoted by f/^(r, tq), is the solution of the equation 

z£^f/:(r,ro) = /i,(r)t/:(r,ro) (8) 

with t/^(ro,ro) = id. 
Notice that 

n(r)f/,^(r,ro) = f/:(r,ro)nM, 

that is, the adiabatic evolution preserves the occupation probability of the band S. 

Now we can adapt to our setting the strong version of the quantum adiabatic theorem, as 
stated in ll27ll . 

Theorem 3.1: Assume that H{u) = H0+U1H1+U2H2 satisfies (HO), and that S is a separated 
discrete spectrum on u; C M^. Let / = [tQ,tf], u : I tu he a curve and set h{t) = H(u{t)). 
Then G C^(J, C-ili)) and there exists a constant C > such that for all r,ro E I 

\\U%T, To) - W,{t, To) II < (1 + |r - rol) . (9) 

Remark 3.2: If there are more than two parts of the spectrum which are separated by a gap, 
then it is possible to generalize the adiabatic Hamiltonian in the following way ( ETI ): 

hair) = h{T) - ieJ2Po.{T)Pa{r) 

a 

where each Pair) is the spectral projection associated with a separated portion of the spectrum, 
partitioning it as a varies. 

Remark 3.3: In general the adiabatic theorem is stated for a time dependent Hamiltonian h{t) 
satisfying the following hypotheses: it is assumed that all the Hamiltonians h{t) have a common 
dense domain V and that the function t H- h{t) is C^{I) and bounded as a function from / to 
Csa{T^,y-), where denotes the space of bounded self-adjoint linear operators from V 

to Ti and V is endowed with the norm of the graph of h{t), for some i E I. These hypotheses 
are satisfied for an Hamiltonian of the form h(t) = Hq + ui(t)Hi + U2{t)H2 under assumption 
(HO), provided that the curve u(-) is C^. 
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In this paper we are particularly interested in the adiabatic evolution across conical intersec- 
tions between eigenvalues. A result in this direction can be found in ETl Corollary 2.5]. In the 
language of control theory it reads as follows. 

Proposition 3.4: Let S : u {Ao(u), . . . , Afc(u)} be a separated discrete spectrum on to. Let 
u°, u^, Uj E cu, j = 0, . . . , k — 1. Assume that Xi{u^) and Ai(u^) are simple for all / = 0, . . . , /c, 
and that, for any j = 0, . . . , /c — 1, Uj is a conical intersection between Xj and Xj+i, with Xi{uj) 
simple if / ^ j,j + 1. Let 7(-) : [0, 1] — )■ be a curve with 7(0) = u° and 7(1) = u-^ 
and such that the eigenstates corresponding to the eigenvalues A; can be chosen along 7 
for all / = 0, . . . ,k. Assume moreover that there exist times < to < ■ ■ • < ^fc-i < 1 with 
7(t,) = Uj, 7(t,) 7^ 0, j = 0, . . . , /c — 1, and that for any / = 1, . . . , /c Xi{'y{t)) is simple for 
every t ^tj, j = 1, . . . , A; — L 

Then there exists C > such that, for any e > 

||7/;(l/5)-e^%|| <Cv^, (10) 

where ^9 G M and ?/'(■) is the solution of equation ^ with -0(0) = 0o(u°) corresponding to the 
control u : [0, — a; defined by u(t) = 'y{et). 

In this paper we are interested in finding control paths along which we have a knowledge 
of adiabatic evolution finer than in (fTOl ). This allows also richer control strategies than those 
described in Proposition 13.41 as it is needed to prove spread controllability. For this purpose we 
write an effective Hamiltonian describing the dynamics inside a two-dimensional band, possibly 
with conical intersections. 

Let us then consider the band constituted by the eigenvalues Aj,Aj+i G S; we can find an 
open domain u' C uj such that {Aj, Aj+i} is a separated discrete spectrum on uj' . 

As above, we consider a control function u(-) G C^(/, w'), for a given time interval /. We can 
then apply the adiabatic theorem to the separated discrete spectrum S' : u i-!- {Aj(u), Aj+i(u)}, 
u G uj': we call P(r) the spectral projection on the band {Aj(u(r)), Aj+i(u(r))} and ]HI(r) = 
P(r)'H its range, which is the direct sum of the eigenspaces of Aj(u(r)) and Aj+i(u(r)). We 
consider the adiabatic Hamiltonian /?.a(r) = /i(r) — z£:P(r)P(r)— i£:P-'-(r)P-'-(r) and its associated 
propagator f/^(r,ro). 

We are interested in describing the dynamics inside EI(r). Since ]HI(r) is two-dimensional for 
any r, it is possible to map it isomorphically on and identify an effective Hamiltonian whose 
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evolution is a representation of f/a(T, ro)|H(ro) on C^. 

Let us assume that there exists an eigenstate basis {</)Q,(r), 0^(7")} of H(r) such that (pa{-), 
belong to C^{I,'H). We construct the time-dependent unitary operator W(r) : IH[(r) — )■ by 
defining for any ip E IHI(r) 

U{T)ij = ei(0,(r),z^) +e2(0/3(r),^), (11) 

where {61,62} is the canonical basis of C^. 
We then define the effective propagator 

U!s{r,To) =U{r)U',{r,To)W{To). (12) 
It is easy to see that f/^g(r, tq) satisfies the equation 

^e^U^Ar, To) = H^ArW^Ar, To), U^siro, Tq) = id, (13) 
where HI^{t) is the effective Hamiltonian whose form is 

HI^{t) = U{T)ha{T)W{r) + ieU{T)U*{T) 

^(Kir) \ _^J{Mr).Mr)) (0^(r),4(r)) 
V X^{r)) '\{Ur),Mr)) (Mr), Mr)) 
Theorem 13.11 implies the following. 
Theorem 3.5: Assume that {Aj,Aj+i} is a separated discrete spectrum on cu' and let u : 
[to,tf] — 7- be a curve such that there exists a C^-varying basis of M{-) made of eigenstates 
of h{-). Then there exists a constant C such that 

II (t/^(r,ro) - W*(r)?7,^ff(r,ro)W(ro))P(ro)|| < Ce{l + \t - to\) 

for every r, Tq G [to,^/]- 

B. Regularity of eigenstates 

Classical results (see [25|) say that the map u 1— P^, where is the spectral projection 
relative to a separated discrete spectrum, is analytic on cu. In particular, eigenstates relative to 
simple eigenvalues can be chosen analytic with respect to u. 

Similar results hold also for intersecting eigenvalues, provided that the Hamiltonian depends 
on one parameter and is analytic. In particular, if S is a separated discrete spectrum on cu and 
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u : J — 7- (jj is analytic, then it is possible to find two families of analytic functions : / — R and 
$j : / "H, j = 0, . . . ,k, such that for any t in / we have S(u(t)) = {Ao(t), . . . , Ak{t)} and 
($o(^)) • • • ) ^kit)) is an orthonomial basis of corresponding eigenstates (see [fT6l . [|25l Theorem 
XIL13]). 

Moreover, we can easily find conditions on the derivatives of the functions Ai,^i: indeed, 
consider a curve u : / — )• such that there exist two families of functions : / — )• M 
and : I ^ v., I = 0, . . . ,k, which for any t E I, correspond to the eigenvalues and the 
(orthonormal) eigenstates of H(u(t)). 

By direct computations we obtain that for alH G / the following equations hold: 

Ai{t) = {Mt)H, + U2{t)H^) ^i{t)) (15) 

(A™(t) - Ai{t)) ^m{t)) = {ui{t)Hi + U2{t)H2) (16) 

An immediate consequence of (fTS]) is that the eigenvalues A; are Lipschitz with respect to t. 
Let u be a conical intersection between Aj(u) and Aj+i(u). Consider the straight line reit) = 
u + t{cos9, sin 6*), t > 0. Then (fT6l) implies that 

lim {(j)j{re{t)), {cosOH, + sin ^i/2)0i+i(re(t))) = 0. (17) 

IV. Conical intersections 

From now on, we assume that the Hamiltonian satisfies hypothesis (HI). Following Re- 
mark 12.11 we always choose the eigenfunctions of H{u.) whose components are real with 
respect to the basis {xi}i defined in hypothesis (HI). In particular, this ensures that the values 
(0i(u),ifo0m(u)), (0i(u),ifi0^(u)) and (0z(u), if20m(u)), /,m = 0,...,/c, are real for any 
u. 

In this section, we investigate the features of conical intersections and provide also a criterion 
for checking if an intersection between two eigenvalues is conical. First of all we notice that 
Definition 12.41 can be reformulated by saying that an intersection u between the eigenvalues \j 
and Aj+i is conical if and only if there exists c > such that, for every straight line r{t) with 
r(0) = u, it holds 

d r 1 

> c. 



^ A,+i(r(t))-A,(r(t)) 



dt 



t=o+ 



Moreover, the following result guarantees that (|6j holds true in a neighborhood of a conical 
intersection. It follows directly from the Lipschitz continuity of the eigenvalues. 
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Lemma 4.1: Let u a conical intersection between Xj and A^+i. Then there exists a suitably 
small neighborhood [/ of u and C > such that 



Aj+i(u) - Xj{u) > C\u - u|, Vu G f/. 



(18) 



Let us now introduce the following matrix, which plays a crucial role in our controllability 
result. 

Definition 4.2: Let ^i,4'2 be a pair of elements of Ti. The conicity matrix associated with 

(^1,^2) is 

' (V^l,i/l^/^2) |((^2,i^l^2)-(^l,i^lV'l)) 

Lemma 4.3: The function {'ipi,ip2) ^ \det Ai{ilJi,ip2)\ is invariant under orthogonal trans- 
formation of the argument, that is if {il)i,il)2Y = (^{'^i,'^2)'^ for a pair tpi,tlJ2 of orthonormal 
elements of H and O G 0(2), then one has | det Miipi, 1^2)1 = I det M {1^1,1/^2) \- 



M{i)i,ip2 



(19) 



Proof: We set = ( „^^°Jina /ros"a )' where ^ = ±L A direct computation shows that 
A4(V'l,^2) = M{ilJi,iJ2) 



COS 2a 
sin 2a 




which immediately leads to the thesis. 



□ 



The following result characterizes conical intersections in terms of the conicity matrix. 
Proposition 4.4: Assume that {Aj, A^+i} is a separated discrete spectrum, and Aj(u) = 
Aj+i(ii). Let {^1,^2} be an orthonormal basis of the eigenspace associated with the double 
eigenvalue. Then u is a conical intersection if and only if A^(?/^i, ?/^2) is nonsingular. 

Proof: Let ro{t) = u + t(cos ^, sin^) and let 0^,0^^^^ be the limits of (f)j{re{t)), (j)j^i{ro{t)) 
as t — 7- 0^ (recall that the eigenfunctions (pj,(pj^i can be chosen analytic along for t > 0). 
Assume that for any e > there exists such that 

d 



dt 



t=0+ L 



that is, by dH, cos {{(Pf,H,(Pf) - H,(P%,))+sm9, {{(Pf,H24>f) - H2<P%,)) < 

e. Moreover, by ([l7]), we have that cos9e{(l)f , Hicpf^^) + sm9e{(f)f , H2(l)%i) = 0. Since 



detA^( 



det 



cos 6p sin 6^ 



sin 9 s cos 9 s 



Mi 



"3 ' '^i+i^ 



<2£(||i/i| 



2IIJ, 
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then, by Lemma |43] and the arbitrariness of e, we get that A^(z/;i,^2) is singular. Thus u is a 
conical intersection when M.{ipi,il)2) is nonsingular. 

Let us now prove the converse statement: assume that u is a conical intersection and, by 
contradiction, that A^(0^,0^^^) is singular, where (pp(f)^j_^_i are defined as above. By definition 
of conical intersection, we have 

cos (3 ((0,^ H^<P^) - + sin (3 ((0^^ H^cj^^) - {ct>U^ H,ct>U)) ^ 0, (20) 

for every /3 e M. 

By ([ni) and (l2Ql) with (3 = 6,ii turns out that the two columns of the matrix A^(0^,0^_^i) 
are not proportional. Thus A^(0^, 0^_,_i) can be singular only if its first column is null. 

For any angle (3, there exists an orthonormal matrix O = ( ^"f^a ?os2 ) such that (0^, (t>j+iV = 
O(0^,0^+i)^ and, calling W = cos/3ifi + sin/Jifa, we have (by (flTl) ) 

= (0^^, iy0^+i) = (cosa^ - sina2)(0^^, iy0,Vi) + sinacosa ((0^%i, 1^0,'+i) - (0,', Wcp'^)) = 
= sin a cos a ((0^^+i, W<j)'^^,) - {<P% 1^0^^)) . 

If (0^^_i, Vr0^+i) - (0^, 1^0^) = 0, the matrix ((0f , 1^0^))z,m=jj+i is diagonal and proportional 
to the identity. Hence the same is true for {{(p^ ,W(j)^))i^rn=j,j+i- This contradicts (|20] ). so that 
it must be sin a cos a = 0, that is, the limit basis is unique and therefore it must be equal to 
{0^,05+1} (up to phases). 

Let us now consider the straight line rp with 

(0,V,gi0,V)-(0^^gi0p 

'^^^ (0^1,i^20,V)-(0,^i/20p- 

Since, as proved above, the limit basis along rp is {0^,0^_,^i} we have that (0^, (cos /^iiTi + 
sin/3iy2)0j) = (0j+i, (cos/3ifi + sin /3if2)0^+i)- By ([B]), this contradicts dJO]), proving that 
A^(0^,0^^]^) is nonsingular. □ 

As noticed above, for any analytic curve that reaches a conical intersection it is possible to 
choose continuously the eigenstates along the curve. A peculiarity of conical intersections is 
that, when approaching the singularity from different directions, the eigenstates corresponding 
to the intersecting eigenvalues have different limits, and the dependence of such limits from the 
direction can be explicitly computed, as shown in the following result. 



January 12, 2013 



DRAFT 



16 



Proposition 4.5: Let u be a conical intersection between the eigenvalues Aj,Aj+i and let 

0°, 0°+! be the limits as t — )■ 0+ of the eigenstates (j)j{ro{t)), </)j+i(ro(t)), for ro(t) = u + (t, 0). 
Consider, for any a G [0,2%), the straight line ra{t) = u + (t cos a, t sin a). Then, up to a sign, 
the eigenstates (f)j{ra{t)), (j)j^i{ra{t)) have limits 

0° = cos H(a)0O + sin H(a)05+, (21) 
= - sin S(a)</)° + cos S(a)0°_^i, (22) 

where H is a monotone function defined on [0, 2n) with S(0) = 0. Depending on the initial 
choice of 0°, (p^^i the range of S is either [0, tt) or (— tt, 0]. Moreover, H(-) satisfies the equation 



cos a, sma 



cos 23(a) 
sm2'E(a) 



-^(0;>0;+i) I 1=0. (23) 



Proof: Let us write 0", 0°^_]^ as in (|2TI)-(|221). Then 3(a) satisfies 

= (0J, (cosaifi + sinaif2)0jVi) 
= cos 23(a) (0°, (cos aHi + sinaif2)0j+i) + 

+ i sin 23(a) ((0j+i, (cosaiiTi + sinaiir2)0j+i) - (0°, (cosai^i + sinai?2)0°)) 

cos 23(a 



cosa,sma 



sin 23(a) 

proving (|23]) . Equation (|23T ) has exactly four solutions for any value of a, differing one from 
the other by multiples of 7r/2. By the Implicit Function Theorem, it turns out that each of them 
is a monotone function defined on [0, 2tx). 

We define 3(-) as the one that satisfies 3(0) = 0. We are left to prove that the range of 3 
is [0, tt) or (— 7r,0]. We first observe that when a = tt the possible solutions of equation (l23l) 
are multiples of 7r/2. If |3(7r)| > 7r/2, then by continuity there should exist a G (0,7r) with 
|3(a)| = 7r/2. This is impossible because of equation (l23l) . Thus 3 maps [0, tt] into [0,7r/2] or 
[— 7r/2,0] and, by symmetry, the claim is proved. □ 

Remark 4.6: From Proposition 14.51 it is straightforward to see that it is not possible to define 
continuously the eigenstates 0j,0j_|.i of H{\i) on a closed path that encloses the singularity: 
after a complete turn, a change of sign appears. 
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V. Non-mixing curves 

Throughout this section we will assume that {Aj,Aj+i} is a separated discrete spectrum on 
some open domain to, and that G w is a conical intersection between the eigenvalues. Without 
loss of generality, in the following we always assume that is the only intersection between \j, 

Aj+i in oj. 

Following Section IIII-Al the effective Hamiltonian H^^, defined as in (fT4)) . (approximately) 
describes the dynamics inside the eigenspaces associated with Aj, A^+i, for u slowly varying 
inside cu. 

When integrating the effective Hamiltonian, the second term in (O gives a total contribution 
that a priori is of order 0(1). In particular the contribution of the non-diagonal terms of H^^ 
induce a (a priori) non-negligible probability transfer between the two levels. 

To tackle this issue we consider trajectories satisfying the following dynamical system 

Ul = -{(pj,H2^j+l) ^^^^ 

Notice that the right-hand side of (|24] ) can be taken real-valued under hypothesis (HI). It is 
defined up to a sign, because of the freedom in the choice of the sign of the eigenstates. 
Nevertheless, locally around points where \j ^ Aj+i, it is possible to choose the sign in 
such a way that the right-hand side of (|24|) is smooth, and, from equation (fT6l) . we see that 
(0j(7(t)), 0j+i(7(t))) = along any integral curve 7 of (l24l) . Here and in the following we use 
the notation 0(7(-)) to denote ^(0(7(-))). 

Let now Ti^ be the real Hilbert space generated by the basis {xj}j defined in Remark 12. 1[ 
and let Gr2('H'^) be the 2-Grassmannian of H^, i.e. the set of all two-dimensional subspaces 
of Ti*. This set has a natural structure of metric space defined by the distance d{Wi,W2) = 
\\Pwi — Pw2h where Pwi,Pw2 the orthogonal projections on the two-dimensional subspaces 
Wi, W2. Lemma |43] allows us to define the function 

F : Gt2{H^) M (25) 

W I det A^(t;i, V2)\, 

where {vi, V2} is any orthonormal basis of G Gt2{'H^). It is straightforward to see that F is 
continuous. 
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Let Pu be the spectral projection associated with the pair {Aj(u), Aj+i(u)}. We know from 
Section ITlI-B I that is analytic on to. Therefore u i-t- P^'Hn'H^ is continuous in Gr2{'H^). Let 
now F(u) := \detM{(j)j{u),(j)j+i{u))\. Since F(u) = F{P^nnn^) and by Proposition ED 
we get the following result. 

Lemma 5.1: The function u i— t- F{u) is well defined and continuous in to. In particular F is 
different from in a neighborhood of u = 0. 

Without loss of generality, we assume from now on that F is different from zero on to. 

Lemma 5.2: There exists a C°° choice of the right-hand side of (|24l) in a; \ {0} such that, if 
u(-) is a corresponding solution, then 

d 



dt I 

on a;\ {0}. 

Proof: Observe that 

d_ 

dt 



A,-,i(u(t))-A,(u(t)) =-F(u(i)) (26) 



Aj+i(u(t))-Aj(u(t)) = ui{{(pj+i, Hi(f)j+i)-{(f)j, Hi(f)j))+U2{{(pj+i, H2(pj+i)-{(l)j, H2(pj)). 



This expression, evaluated along the solutions of (l24l) . is equal either to F(u(t)) or to — F(u(t)), 
depending on the choice of the sign in (l24l) . Since F{u) 7^ on w, there exists a unique choice of 
this sign such that equation (|26|) is satisfied. The local smoothness of the eigenfunctions ensures 
that this choice is smooth. □ 

We now define the non-mixing field, denoted by Xp, as the smooth vector field on cu \ {0} 
identified by the preceding lemma. Its integral curves are C°° in uj\ {0}. Moreover, its norm is 
equal to the norm of the first row of J^{(pj, <pj+i), and therefore bounded both from above and 
from below by positive constants in a; \ {0}. 

By considering Aj+i(u) — Aj(u) as a local Lyapunov function, the above results lead to the 
following proposition. 

Proposition 5.3: There exists a punctured neighborhood f/ of such that all the integral 
curves of Xp starting from U reach the origin in finite time. 

Our purpose now is to prove that each of these curves admits a C°° extension up to the 
singularity. As a preliminary result we get the following. 

Proposition 5.4: Let u = be a conical intersection with Aj(0) = Aj+i(0), and let the map 
u I-)- {Aj(u), Aj+i(u)} be a separated discrete spectrum on a neighborhood of 0. Then, for any 
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Ci > 0, there exist a neighborhood / of t = and C2 > such that for any trajectory 7(-) 
with 7(0) = 0, |7(0)| = 1 and ||7lU-(/) < Ci, one has \\M-f{t))\\ < C2, 1= j,j + 1, for every 
tel\{0}. 

Proof: Let us assume without loss of generality that Aj(0) = Aj+i(0) = 0. For t ^ define 
p(t) = 7(t) - 7(t)/t. Notice that p{t) = l^iO)t + o{t). 
By (fT6l) . we have 



-j, {llHl + 72i/2)0j+l) , (Pl^^l + P2^2)0j+l) 



(0.(7(t)), 0,>i(7(t))) = --^^-^ + ---^'--7- .^-^.^.^v . (27) 



Notice that (0j(u), (wiifi + M2i^2)0i+i(u)) = -(0j(u), ifo0i+i(u)) = -(0j(u)-Po0j(u), i^o^vj 
Po0i+i(u))). Since 

||ifo(0i(u) - ^O0j(u))|| = \\\j{u)(j)j{\l) - UiHi(j)j{u) - U2H2(j)j{u)\\ 

< \X,{u)\ + \u\i\m\ + \\H2\\) 

< snp\{(P,{v),{mH^ + U2H2)<l),{v)) I + |u|(||i7i|| + ||i72||) 

< 2i\m\ + \\H2\\)\u\ 

and by smoothness of the projector, we get that | (0j(u) — Po0i(u), Hq (0j_|_i(u) — Po0j+i(u)) ) | < 
2C{\\Hi\\ + ||i/2||)|u|2, for a suitable C > 0. Being |7(t)| = 0{t) and Xj+i{u) - Xj{u) > c\u\ 
(Lemma 14.11) . we deduce that the modulus of the first term in the right-hand side of (ITtT ) is 
uniformly bounded. The uniform bound of the second term is a trivial consequence of the 
fact that |p(t)| = 0{t) and that |7(t)| > c\t\, for some c > 0, if t is small enough. Thus 
I (0j(7(t)), 0j-i_i(7(t))) I is uniformly bounded. 
Let us write P^ = id — Pu. Since commutes with H(u), one has 

(/7(7(t)) - A,+i(7(t))id)P4)0,+i(7(t)) = -P4)(7i(t)//i + 72(t)i^2)0,+i(7W). 

Since H(u) — Aj+i(u)id is invertible on P^H with uniformly bounded inverse on to, we get that 
||P^^)0-,+i|| is uniformly bounded on /\ {0}. Thus we obtain that ||0j+i|| is bounded, uniformly 
on the set of curves 7(-) satisfying the assumptions of the proposition. The same holds for 

□ 

Corollary 5.5: Let u = be a conical intersection with Aj (0) = Aj+i(0), and let the map 
u I— {Aj(u), Aj+i(u)} be a separated discrete spectrum on a neighborhood of 0. Denote by 
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(t>i{p,9) the eigenstate 0/(pcos^,psin^), / = j,j + 1, where {p,0) are angular coordinates 
around 0, i.e. p = \u\ and 9 = arctan^. Set 4)i{0,9) = \imp_^o+ (j)i{p,9). Then the function 
(piip, 9) is continuous in [0, R\ x [0, 27r], for some i? > and / = j, j + 1. 

Proof: If p > 0, the function (f)i(p,9) can be defined continuously. Moreover, the function 
9 I-)- (pi{0,9) is uniformly continuous, thanks to Proposition 14.51 

Let us now consider a sequence {pk,9k) converging to (0,^). Then we have 

\Mpk,Ok)-Mo,o)\ < \Mpk,Ok)-Mo,dk)\ + |0;(o,^fc) -0;(o,^)l 

where Ci comes from Proposition 15.41 and the second term goes to zero as k goes to infinity. □ 



Proposition 5.6: The eigenstates (f)j,(j)j+i can be extended continuously to the singularity 
along the integral curves of Xp, and, in a small enough punctured neighborhood of u = 0, the 
integral curves of Xp admit a extension up to the singularity included. 

Proof: We prove that the scalar product Xp ■ (— ^2, ni)-^/|u| goes to as |u| — 0, that is, 
the tangent to the curve has limit when u approaches zero. This, together with Corollary 15. 5J 
implies that the eigenstates 0^,0^+1 are continuous along the integral curves of Xp, and then 
the vector field Xp itself is continuous along its integral curves, up to the singularity included. 
Therefore its integral curves admit a extension up to the singularity. 

To prove that Xp ■ {—U2, ni)^/|u| goes to as |u| — )■ 0, we show that there exists a constant 
C > such that 

k{u) ■= \Xp{u) ■ (-M2,Mi)| < C\u\\ (28) 

Since k(u) = |(0j(u), {uiHi+U2H2)(j)j+i{u.))\, the thesis comes from the estimates in the proof 
of Proposition 15.41 □ 

We recall that, since integral curves of the non-mixing field Xp are C\ then the spectral 
projection associated with the pair {Aj(u), Aj+i(u)} is along each of them. This permits 
to prove the following result. 

Proposition 5.7: For any integral curve 7 : [—77, 0] — t- u; of Xp with 7(0) = there exists 
a choice of an orthonormal basis of the eigenspace associated with the double eigenvalue 
Aj(7(0)) = Aj+i(7(0)) that makes the eigenstates 0j(7(t)), 0j+i(7(t)) on [—ri,0]. 
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A'p(7(0)) 



Proof: We notice that on the integral curves of Xp the eigenstates relative to the eigenvalues 
Aj,Aj+i satisfy the equation P^i^t)4>j{'^{t)) = P^(t)0j+i(7(t)) = 0, which implies 

P.itMlit)) = 0,(7W) P,(i)0,+i(7(t)) = 0,+i(7(t)) (29) 

for t G [— ?7,0). The thesis follows from the continuity of P^i^t)-,4>j{'l{^))-,4>j+i{'l{t)) on [~'75 0]. 
□ 

Corollary 5.8: Let 7 : [—77, 0] — a; be an integral curve of Xp with 7(0) = 0. Then 7(-) and 
the eigenstates </)j(7(-)), </)j+i(7(-)) defined in Proposition 15.71 are on [—77,0]. 
Proof: Extend Xp{'^{t)) by setting 

-(0,(7(0)), ^r20,+i(7(O))) 
(0,(7(0)), i^i0,+i(7(O))) 
where 0,(7(0)), 0,+i(7(O)) denote the limits of the eigenstates as defined in Proposition 15.71 
Then Xp{^{-)) is on [—77, 0], which implies that 7(-) is on [—77, 0]. We differentiate equation 
(|29l ) to prove that 0,(7(-)), 0,+i(7(-)) are on [—77, 0]. Repeating recursively the argument we 
prove the thesis. □ 
We stress that, thanks to Proposition 15.71 if we define the adiabatic Hamiltonian ha{T) = 
H{-^{t)) — ieP^(r)Pj{T) — «e-P^^)-P^T-), T = et, along integral curves of Xp, then it is possible 
to define the associated effective Hamiltonian, as in equation (fT4l) . 
The following result is crucial to our controllability strategy. 

Proposition 5.9: For every unit vector w in there exists an integral curve 7 : [—77, 0] — )■ w 
of Xp with 7(0) = such that 

at) 

lim - — -— -7 = w. 



Proof: Equation (|24|) rewrites as 

p= (0,(p,^),(-cos^i72 + sineifi)0,+i(p,^)) (30) 

e = ^(0,(p,^),(cos^ifi + sin^i72)0i+i(p,^)). (31) 

On a neighborhood C of the singularity, there exist two constants < ci < C2 such that 
ci < IpI < C2, and the right-hand side of (IBTT ) is bounded from above, by (|28l) . We choose the 
sign of the functions 0,, 0,+i in such a way that p < 0. 

Fix 9 E [0,27r] such that w = (cos^, sin^). Consider, for k large enough, the solutions 
(Pfc(")) ^fc(-)) of (I30l)-(|3T|) with Pa;(0) = 1/k and 6*^(0) = 9, for t belonging to some common 
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interval [—77,0], where 77 > is small enough, in order to guarantee that the solutions do not 
exit from U. By Ascoli-Arzela Theorem, up to subsequences, iPk{-),dk{-)) converges uniformly 
on [— ?7,0] to some {p{-),6{-)). 

In particular, for any r G [—77, 0], {pkir), Ok{T~)) converges in U . By the uniform boundedness 
of p, the range of {pk{-), ^fc(')) [~V^ ''"] contained in a compact subset K C f/\{0} for every 
k. Since the vector field is smooth on K, the curves {pki'), (^k{-)) converge uniformly on [—77, r] 
to the solution of (l30l)-(l3T1) with initial condition p(r) = p(r) and 9{t) = 9{t). Therefore for 
t e [-r], t] (p(-), 9{-)) is a solution of dSQj-dll])- Since r is arbitrary, and ^(0) = lim^ 6k{0) = 6, 
p(0) = limfc pfc(O) = 0, we get the thesis. □ 

We conclude this section by proving a result of structural stability of conical intersections. 

Theorem 5.10: Assume that H{u) = Hq + uiHi + U2H2 satisfies (HO)-(Hl) and let u be 
a conical intersection for H{u) between the eigenvalues Xj and Aj+i. Assume moreover that 
u I— 7- {Aj(u), Aj+i(u)} is a separated discrete spectrum in a neighborhood of u. Then for every 
e > there exists 6 > such that, if H{u) = Hq + UiHi + U2H2 satisfies (HO)-(Hl) and 

11^0 - H4 + 11^1 - + \\H2 -H2\\< S, (32) 

then the operator H{u) admits a conical intersection of eigenvalues at u, with |u — u| < e. 

Proof: Continuous dependence of the eigenvalues with respect to perturbations of the 
Hamiltonian ensures that, if 6 is small, then H admits two eigenvalues Xj, A^+i close to Aj, A^+i. 
Moreover {A^, Aj+i} is separated from the rest of the spectrum, locally around u. Fix now e > 
in such a way that the vector field Xp points inside the ball i?(u, e) at every point of its boundary 
(this is possible because of (l28l) ) and F(u) > c > on B{u,e). If 6 is small enough then 
Xj 7^ Aj+i on dB{u,e). Similarly, since the conicity matrix Ai varies continuously with respect 
to Hi,H2, and by continuity of the function F defined in (l25l) . we can take 6 small enough 
such that IdetA^I > c/2 for any perturbed Hamiltonian. This allows us to define, whenever 
Xj 7^ Aj+i, the non-mixing field Xp associated with H and corresponding to the band {Xj, Aj+i}; 
as in Lemma 15. 2[ we choose Xp in such a way that the time derivative of Aj+i — Xj along the 
integral curves of Xp is smaller than —c/2 and Xp is smooth. In addition, by the uniform 
continuity on dB(\i,£) of the eigenfunctions with respect to perturbations of the Hamiltonian, 
if 5 is small enough, then Xp points inside B(\i,e) at every point of dB{\i,E). 
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Fix an Hamiltonian H{-) satisfying (|32l ). Any trajectory 7(-) of Xp starting from B(u,e) re- 
mains inside -B(u, e) in its interval of definition and reaches in final time a point u corresponding 
to a double eigenvalue Aj(u) = Aj+i(u). The conclusion follows from Proposition 14.41 □ 

VI. Proof of Theorem [23] 

Based on Proposition 15. 9[ we consider below trajectories of the following kind: given a conical 
singularity u and a pair of unit vectors wi, W2 G M^, we concatenate the integral curve of Xp 
arriving at u with direction wi and the integral curve of —Xp exiting u with direction W2. 
Even if we are not using this fact in the paper, it turns out that, if wi = W2, then such curve is 
globally C°°. 

Proposition 6.1: Let u = be a conical intersection between the eigenvalues Aj,Aj+i and 

let 0°, (p'j^-j^ be limits as r — )■ 0^ of the eigenstates 0j(r(r)), </)j+i(r(r)), respectively, for r(r) = 

(r, 0). Let 7 : [0, 1] — )■ c<j be a piecewise C°° curve such that 7(ro) = for some tq G (0, 1), 

7(r) = A'p(7(r)) in [0,ro] and 7(r) = -Xpi-yir)) in [tq, 1]. Define a_, a+ by 

7(t) 'ii'^) 
lim J — i-—r- = — (cosa_, sina_) , lim - — -— -- = (cosa+, sina+). (33) 

Then there exists C > such that, for any e > 0, 

||^(l/e) -pie^^^0,(7(O)) -P2e^''^'+V,+i(7(0))|| < Ce (34) 

where 'dj+i G M, 'tp(-) is the solution of equation ^ with tp{0) = 0j(7(O)) corresponding 
to the control u : [0, — > w defined by u(t) = 'j^et), 

pi = I cos (H(a+) - H(a_)) |, p2 = | sin (S(a+) - H(a_)) |, 



and S(-) is defined as in Proposition! 

Proof: We consider the Hamiltonian H{u{t)), t G [0, l/e]. Since the control function u(-) 
is not at the singularity, we cannot directly apply the adiabatic theorem. Instead, we consider 
separately the evolution on the two subintervals (in time t) [0,tq/e] and [To/e^l/e]. 

Since the eigenstates 0j(u(t)), 0j+i(u(t)) are piecewise we can apply Theorem |331 in order 
to study the evolution inside the space Pu{tyH. We can then construct the effective Hamiltonian, 
which is diagonal on both intervals (in time r) [0, tq] and [tq, 1]. Remark that the operator- valued 
function W(-), defined in equation (fTTI) . has a discontinuity at tq but has continuous extensions 
on both intervals [0, tq] and [tq, 1]. 
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Let (pf = lim^ ± (j)A'^{T)). Integrating the effective Hamiltonian we get 



f/,^(ro, 0)^(0) = e^V- 
for some 9? G M . By Proposition 14.51 we have 

0+ = COS {^{a+) - 0" + sin - 07+1 

0^+1 = - sin {'&{a+) - '&{a^)) (p] + cos {'&{a+) - i^ia-)) </)J+i. 
Then, since the effective Hamiltonian is diagonal, we get 

[7,^(1, 0)m = e*''^ cos (^(a+) - ^(a_)) 0,(7(0)) + 6*^^^+^ sin (^(a+) - ^(a_)) 0,+i(7(O)), 
and then, applying the adiabatic theorem, 

ll^(r) -pie^'^^- 0,(7(0)) -p2e^'^+V,+i(7(0))|| < 

where C is a constant depending on the gap and on 7. □ 
Remark 6.2: For control purposes, it is interesting to consider the case in which the initial 
probability is concentrated in the first level, the final occupation probabilities pi and are 
prescribed, and there is an integral curve of Xp connecting u° to the singularity. Except for the 
special cases pi = 0,1, there are exactly two integral curves of —Xp starting from the singularity 
that realize the required splitting (in the sense of Proposition 16.11) . 

Choosing /3 G [0,7r/2] such that (pi,P2) = (cos /3, sin /3), we obtain that the two possible 
values for a+ are 

a+ = (/3 + H(a_) + fc+vr) a+ = H"^ (-/3 + H(a„) + fc^vr) , 

where k^, E Z are chosen in such a way that /3 + H(a;_) + /c+tt and —(3 + H(a_) + k^Tv 
belong to the range of 3. 

If {pI,pI) = (0, 1), then the path is unique with a+ = a„ + vr, while if {pi, pi) = (1, 0), then 
the unique path satisfies a+ = 

Proof of Theorem I2.5i For simplicity, we consider the case in which ^(0) = 0o(u''). The 
general case can be treated similarly. 

Recall that for any conical intersection between two eigenvalues of a separated discrete 
spectrum there exists a neighborhood of the intersection where the two eigenvalues are well 
separated from the rest of the spectrum. Let us consider these neighborhoods for the intersections 
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xj, j = 0, . . . , /c — 1, and let us call them uj. Define on each uj \ {xj} the vector field as 
in Section |Vl 

We construct the path 7(-) as described below. 

First choose a smooth path cro(-) starting from u° and reaching cuq along which all the eigenval- 
ues in S are simple. Concatenate ao with an integral curve of Xp that reaches the point xq- Then 
choose a° as one of the angles realizing, for the two-levels system associated with the energy 
levels Aq, Ai, the splitting from (1, 0) to (p^, 1 — p^), as explained in Remark and continue 
the path with the integral curve of —Xp with outgoing tangent parallel to (cosal, sina° ). 

Join the latter with a smooth path (Ji(-) connecting loq to cui along which all the eigenvalues 
in S are simple, and then prolong it with an integral curve of Xp that reaches the point Xi. As 
above, compute an angle a\. that realizes the splitting (for the two-levels system associated with 
the energy levels Ai, A2) from (1 — pj, 0) to {pi, 1 — {pf +pl)), and, as above, continue the path 
with the integral curve of —Xp with outgoing tangent parallel to (cos a^, sin a^). 

Repeat this procedure iteratively until the required spread is realized. Then reach the final 
point with a path along which all the eigenvalues are simple. We assume without loss of 
generality that the final time is equal to one. 

For e > 0, consider the Hamiltonian H{u{t)) = H{'j{et)), and set r = et. 

As long as 7(r) G \ U^lQUi, we approximate the dynamics of H{u) using the adiabatic 
Hamiltonian 

k 

K{t) = H{^{T))-isY,Pi{r)Pi{r)-iePi{T)Pi{T) (35) 

1=0 

where Pi{t) is the spectral projector onto the eigenspace relative to \i{^{t)) and -Ps"('r) = 
id-Ef=o^K^). 

The evolution associated with (l35l) conserves the occupation probabilities relative to each 
energy level in S, therefore the evolution of iJ(7(r)) approximately conserves these occupa- 
tion probabilities, with an error of the order e, as prescribed by the adiabatic theorem (see 
Remark Il2l) . 

For 7(r) G Uj, j = 0, . . . , — 1, we use instead the adiabatic Hamiltonian 

k 

K{r) = H{^{t)) - zeP,.,+i(r)P,,,.,i(r) - Z5 ^ Pi{r)Pi{r) - ieP^{T)P^{T) (36) 

(=0 

where Pjj+i(r) is the spectral projector relative to {Aj(7(r)), Aj+i(7(r))}. 
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The evolution associated with (|36l) conserves the occupation probabilities relative to the band 
{Aj, Aj+i}, to any other energy level in S and to its remainder in the spectrum. Moreover, thanks 
to the choice of the field Xp, we can also compute the evolution given by (|36l ) inside the band 
Xj+i} (which is the one described in Proposition 16. II ). 

We end up with final state satisfying 

k 

\m/e)-J2pie''''U^')\\<Ce, 

1=0 

for some i^q, . . . ,'dk E M and some C > determined by the adiabatic approximation. Thus the 
system is approximately spread controllable and the theorem is proved. □ 

VIL Mildly mixing curves 

In the previous section we constructed some special curves along which the effective Hamil- 
tonian has a simple form, whose evolution is quite easy to predict. In this section, we consider 
more general curves passing through the singularities. 

We prove below a variation of Proposition 16.11 which generalizes to broken curves the result 
in [|27l Corollary 2.5]: if we choose any piecewise regular curve with a vertex at the conical 
singularity, then we obtain a distribution of probability between the two levels similar to the one 
described by Proposition 16.11 In this case, if the final time is 1 /e, the error is of order ^/£. 

Moreover, we prove that the integral curves of Xp are not the only ones that realize the best 
accuracy (that is, an error which is of order e for a final time equal to l/e): indeed, this can 
be obtained with any curve whose first and second derivatives at the singularity are the same as 
those of an integral curve of Xp. 

Let us consider a curve 7 : [0, tq] — t- lo such that 7(ro) = corresponds to a conical 
intersection between Xj and A^+i, and ^(tq) ^ 0. Assume moreover that (f)j,(pj^i are along 
7 (recall that this is true for analytic curves). Let us consider the Hamiltonian H{'j{et)),t G 
[0,ro/£:], and the adiabatic Hamiltonian (l36l) . Up to a factorization of the trace, the effective 
Hamiltonian reads 

a(r) —ieb{r] 
iehir) —aij) 



where 



_ V.(7(r))-A,(7(r)) ^ ^^^^ ^ {*,.l(7(r))>,(7(r))), 
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and the dynamical system associated with H^g is 

zez = Hl^z (37) 

where z = 1 (as usual, here and below the dot indicates the derivative with respect to 

r). 

Note that the condition of conical intersection implies the existence of two positive constants 
Ci, C2 such that, for r close to tq, —C2 < ^(t) < — Ci (as a consequence Ci\t — tq\ < a{T) < 
C2\t — To\). As for &(r), it is by hypothesis. 

We set D{t) = ^"^'^J/^ „a(°)/£ U{t) = exp (— z D{s)ds), and we perform the change of 
variable ( = Uz, so that ( evolves according to the dynamical system i( = H^g{T)(, where 

1 , . 1/^0 -ib(T)eT lo ^i')'^' 

H!,{r) = -U{r)H:,{r)U{rr' + U{r)Uirr=\ ^ 

^ \ ib{T)e~Tfo «(^)'^^ 

Let us express the evolution operator for H^f^ in the form 



M,(r,0) 



iy{T) 11 



We claim that || Me(r, 0) - id|| < C^/e, for some C > 0. 
From z^Me(r, 0) = ^(r)M£(r, 0) we get the equations 

v = fib 

fi = ^fia — ub 

with initial data z/(0) = l,/i(0) = 0. Since det = 1 we have that \u\ and \ii\ are bounded, 
and then, from the boundedness of b, we get that also |z>| is bounded. 

We recall that along the integral curves of Xp the effective Hamiltonian is diagonal and its 
evolution is exactly f/(r), so that the equations above are solved by z/ = 1, /i = 0. If b{-) is 
not identically equal to zero, then the evolution is not exactly diagonal, but it mixes the two 
components of C (and z). The error done by assuming the evolution diagonal can be estimated 
by evaluating the term fi. 

By variation of constants we have 



/i(r) = - r eTl>(-)'i-b{s)u{s)ds. 
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We can rewrite 

/"(t-q) = / ds 

a[s) 







2i\ J a{s) lo Jo 2i\ J ds \ a{s) J 
e¥r««--l £/3;;o ,^ 6(0) 



2z r-i>ro a(r) 2z V / a(0) 



2i L \ J ds\ a(s) 



ds. (38) 



Since d(ro) 7^ 0, we obtain that lim^^^Q ^ " * ^^^^ — — = 0. The second term in the equation 
above is of order e. Then, we are left to estimate the integral term. We have 

Jo V / ds \ a[s) J Jo a^{s) 

We consider the change of variables s ^^(s) = | J^" a{r)dr so that 

^0 |eT/?"{'^)*- - 11 e /-^-W |e*^-l| 

ds = - ^ , ^ , , , , (ix. 



a2(s) 2 Jo a^{^-^{x)) 

From the estimates above on a we easily get that a^{^~^{x)) > Ce^l'^x^l'^ for a suitable positive 
constant C. 

Since Jg^°° -^^^3^6/2; < Jp^°° minj^"^/^, 22;~^/^}(i2 < +00, we immediately obtain that the 
integral in (|38T ) is of order Therefore, ||Me(ro,0) — id|| is of order 

If 7 is defined also for r > tq and is globally C^, we recover Corollary 2.5 in [57]. If, instead, 
7 is continuous and piecewise C^, with different tangent directions at the singularity, then we 
can repeat the same argument as in Proposition 16.11 at the singularity the limit basis rotates 
instantaneously and we consider separately the evolution of two different adiabatic Hamiltonians. 
The rotation of the limit basis spreads the probabilities as described by equations (|2TI) - (|22)) . and 
this leads to a controllability result in the spirit of Theorem 12. 5[ where the error is of order ^ 
if the final time is Xje. 

The following result shows that the value of 6(r) at the instant where the curve attains the 
singularity depends only on the 2-jet of the curve at the singularity. This allows us, using 
piecewise analytic curves that have the same 2-jet at the singularity as an integral curve of Xp, 
to obtain a controllability result equivalent to Theorem 12.51 (see Proposition 17.21) . 
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Lemma 7.1: Let 7(-) and 7(-) be two curves on cu such that 7(ro) = 7(ro) = 0, where 
is a conical intersection between Xj and A^+i, with 7(ro) = 7(70) 7^ and 7(ro) = 7(7-0). Let 
r/(r) = (0j+i(7(r)),0i(7(7-))) and r/(r) = (0j+i(7(r)), 4-(7(r))). Then lim,^,„ (r/(r) -r/(r)) = 
0. 

Proof: First of all, we remark that |7(t)— 7(r)| = o((r— tq)^) and |7(t)— 7(r)| = o{\t—tq\). 
As in the proof of Proposition 15.41 we can prove that there exists C > such that for 
every point u G \ {0} and any unitary vector w G the directional derivative along w 
satisfies ^ C'/I^l- Then we obtain that ||0i(7(''")) ~ = o{\r — tq\), I = 

j,j + 1. Moreover, by (1151) we know that the eigenvalues are Lipschitz in a neighborhood of the 
intersection. 

From (|28T ) we have 

_ _ ((0.(7(r)) - 0,(7(r))) , (71^1 + 72i^2) 0,+i(7(r))) 



+ 



A,+i(7(r)) 


-A,(7(r)) 






(72 - 72)^2)^ 


^.+i(7(r))) 


A,+i(7(r)) 


-A,(7(r)) 






(0,+i(7(r)) - 


</'.+i(7(r)))) 



A,+i(7(r))-A,(7(r)) 

1 1 



+ (0.(7(r)), (7ii^i + 72^^2)0,+i(7(r))) 



,A,+i(7(r)) - A,(7(r)) A,+i(7(r)) - A,(7(r)) 
By previous estimates it follows that all the terms in the right-hand side of the equation above 
go to zero as r goes to tq. □ 

Proposition 7.2: Let u = be a conical intersection between the eigenvalues Aj,Aj+i and 
let 0°, (f)^j^^ be limits as r — )■ 0+ of the eigenstates (f)j{r{T)), 0j+i(r(r)), respectively, for r(r) = 
(r, 0). Let 7 : [0, 1] — (jj be a curve such that there exists tq G (0, 1) with 7(r) = if and 
only if r = tq, 7 analytic on [0,tq] and [tq, 1], and 7(''o'^) 0- Let a_ and a+ be the angles 
describing respectively the inward and the outward tangent direction at the singularity, as in 
(l33l) . Assume that the integral curves of Xp having the same inward and outward tangents as 
7 at the singularity possess also the same 2-jet as 7 at the singularity. Then there exists C > 
such that for any e > 

11^(1/5) -pie^''^ 0,(7(0)) -p2e^''^+V,+i(7(0))|| < Ce (39) 
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where ^9^+1 e M, 'ip{-) is the solution of equation & with 'ip{0) = 0j(7(O)) corresponding 
to the control u : [0, — a; defined by u{t) = 7(£t), and 

Pi = I cos (i9(a;+) — | P2 = \ sin — |, 

with defined as in Proposition 14.51 

Proof: By Lemma ITTI the function 6(r) = (0j+i(7(T)), 0j(7(t))) goes to zero as r goes 
to tq. Moreover, the analyticity of 7(-) easily implies that the term ^(6(s)i^(s)/a(s)) appearing 
in (l38l) is bounded. Thus |/i(ro)| < Ce for a suitable C > 0. 

Then \\Mi;{tq,0) — id|| is of order e, where is the evolution operator defined above. We 
can obtain an analogous estimate for || M£(r, Tq) — id||, r > tq. This completes the proof. □ 
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